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A quantum mechanical picture is presented to describe the 
behavior of confined spinous in a variety oi S = 1/2 chains. 
The confinement is due to dimerization and frustration and 
it manifests itself as a nonlinear potential V{x) oc la;!', cen- 
tered at chain ends {b < 1) or produced by modulation kinks 
(& > 1). The calculation extends to weak or zero frustration 
some previous ideas valid for spinous in strongly frustrated 
spin chains. The local magnetization patterns of the confined 
spinous are calculated. A (minimum) enhancement of the lo- 
cal moments of about 11/3 over a single S — 1/2 is found. 
Estimates for excitation energies and binding lengths are ob- 
tained. 

PACS numbers: 64.70.Kb, 75.10.Jm, 75.50.Ee 



I. INTRODUCTION 

Dimerized and frustrated spin chains have attracted 
considerable interest in recent years. This is due to the 
recent emergence of a variety of experimental quasi-one- 
dimensional systems containing locahzed electrons that 
can be described by spin chains. Among these systems 
there is a sizable number which are gapful due to dimer- 
ization in their low temperature phase, such as CuGeOa 
1 1, a'-NaVsOs §, (VO)2P207 §, Cu(N03)2 • 2.5H2O 
||, CUWO4 §, and Cu2(C2Hi2N2)2Cl4 §. The first 
two dimerize as a consequence of the coupling to the lat- 
tice degrees of freedom, whereas the other compounds 
are intrinsically dimerized, i.e. the dimerization does not 
depend on temperature. The effect of doping on these 
substances represents an interesting issue, in particular 
since it pertains also to defects such as missing spins, 
broken chains and chains of finite length. 

As a general rule, defects in low-dimensional antifer- 
romagnetic spin systems which can be described by a 
RVB-type ground state 0] involve S* = 1/2 states in their 
vicinity The appearance of a certain impurity spin 
at the edges of spin chains can also be discussed in the 
framework of a nonlinear a model for general spin S . 
We restrict ourselves to S* = 1/2, weakly dimerized spin 
chains. The basic idea there is that without the defect 
each spin has a partner with which a singlet is formed. If 
the defect removes one of the partners, the other becomes 
a free S = 1/2 spin which we will henceforth call spinon 



to distinguish it from the other singlet forming spins. We 
understand that this spinon comprises also some dressing 
of the bare 5=1/2 spin, i.e. it is not localized at just one 
site. Since the singlet pairing is not static in the RVB- 
picture the spinon is able to move some distance away 
from the defect. In a critical, gapless system, it will be 
delocalized. This delocalization, however, disappears as 
soon as the couplings are modulated. 

It is the purpose of this work to elucidate in which way 
an explicit dimerization acts as a confining potential for 
the spinon motion leading to its localization. The generic 
Hamiltonian reads 



H = J 



(1 + {-sy)SA+i 



i+2 



(1) 



where S parametrizes the dimerization and a the relative 
frustration by next-nearest neighbor coupling. Previous 
works, e.g. ]To|-[l3t, viewed dimerization already as con- 
fining potential. We understand that there is no confine- 
ment without dimerization = 0) . Here we will develop 
a quantum mechanical picture that treats the cases of 
small and large frustration on equal footing. Although 
this picture is not exact in all details, nevertheless it is 
able to reproduce the main features of the problem on a 
semi-quantitative level. The physical quantities consid- 
ered here are binding energies and local magnetizations. 
In addition, it will be shown that the confinement de- 
pends on the degree of frustration and it will be sublinear 
in the region of low frustration. 

For the sake of concreteness, the case of a chain end is 
illustrated in Fig||. 
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FIG. 1. Three possible configurations at chain ends. The 
thick (thin) soUd lines stand for strong (weak) bonds; the open 
eyelets stand for singlets, the arrow for an unpaired spin, (a) 
weak bond at the chain end corresponds to a free spin; (b) 
strong bond at the chain end corresponds to no free spin; (c) 
configuration after two hops of the free spin from (a). Note 
the misaligned first two singlets at the strong bonds. 

Dimerization localizes spin singlets mostly at the 
strong bonds. If in a dimerized chain one spin is missing 
its singlet partner is freed. If the missing spin had a weak 
bond to the left (assuming without loss of generality hor- 
izontal chains) and a strong bond to the right, the free 
spin situation corresponds to the one in Fig.|l|(a). The 
configuration in Fig.^b) is then the reflected configura- 
tion found on the left of the missing spin. If the missing 
spin had a strong bond to the left and a weak bond to 
the right, the free spin situation is found to the left of 
the missing spin. The configuration in Fig.|l|(a) is then 
a reflected image whereas the configuration in Fig.|l](b) 
is found to the right of the missing spin. Fig.|l|(c) illus- 
trates the situation after two hops of the free spin away 
from its origin in (a). The crucial point is that the sin- 
glets to the left of the spin in Fig.|l|(c) are no longer at 
the strong bonds. This implies an energy loss which can 
be viewed as an attractive potential which ties the free 
spin to its origin. Besides chain ends also solitonic mod- 
ulations, such as kink defects, will be considered in this 
paper. They can also be viewed as spinon traps with 
confining potentials. 

The article is organized in the following way. In the 
next section, the short-range RVB spinon states that de- 
termine the main part of the correlations will be intro- 
duced. Subsequently, the forms of the kinetic and the 
potential energies which govern the spinon dynamics will 
be discussed, i.e. the Schrodinger equation of the prob- 
lem will be setup. In the fourth and the fifth sections 
the quantitative results will be discussed and compared 
against computational calculations for chain ends and 
for kink-modulated chains, respectively. A summary will 
conclude the article. 
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for which we aim to obtain a continuum description. This 
means that we assume that for the low energy behavior 
it is sufficient to treat a.i as a slowly varying function 
of i. Note that we introduced the factor (—1)' to focus 
on the energetically low-lying states. This can be seen 
most easily for a = l/2[|6jl4,l^]. In order to define a 
normalized wave function i]j{x) let us calculate the norm 
oi\v) 




(4) 



The approximate step is valid if varies in fact slowly 
with i. Eq.(^) tells us that passing to a normalized 
continuous wave function ip{x) means i — > x/2 and 
Oi y/2/^->p{x). The factor 2 introduced here enables 
the replacement of x by the actual site numbers at the 
end of the calculation. 

Now we calculate the local magnetizations mi :— (S'f). 
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A contribution is found only if S'f is situated between the 
two spinous which equals one half of the overlap (||) . For 
I odd one finds 



{v\Snv) « |a(/-i)/2|^ 
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II. STATES: NORM AND MAGNETIZATION 

Let us denote by \i) an up-spinon at site 2z -I- 1; the 
other spins are all paired to nearest-neigbor singlets. In 
this convention, the state in Figj^(a) is denoted by |0) 
and the state in Fig.|l|(c) is denoted by |2). These states 
are not orthogonal but their overlap [p^PJlS]] is given by 



(2) 



This overlap arises from a Neel type sequence of up-spins 
and down-spins between the two spinous. Let us consider 
a state 



a 



(;-i)/2l 



(6) 



where n and m stand for half the distance to site / on the 
left and on the right. We do not take the chain end into 
account, i.e. we ignore the fact that summation for small 
I is truncated on the left. For even sites I one obtains 



(v\snv)^-l\ 



0-1/21 
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With the above calculation the local magnetizations have 
been linked to the probability of finding the spinon at a 
given site I. The sum J^i "^i equals to one half as it has 
to be for a global S = 1/2 state. Note that to obtain this 
result one has to sum even and odd sites separately, and 
that I changes, thus, by two from site to site. 

Summing the moduli |m;| one gets immediately 
i . This means that the antiferromagnetic correlations 
induce total local moments that correspond to the mo- 
ments of 11/3 = 3.667 independent spins S = 1/2. This 
enhancement factor of 11/3 illustrates why even a low 
concentration of dopants may induce considerable anti- 
ferromagnetism in a compound |]l7|,^,p^ . Moreover, this 
enhancement is not due to criticality but it is already 
built-in in the short-range RVB state. The possible co- 
existence of dimerization and alternating local magne- 
tization was nicely demonstrated by Fukuyama et al. 
by analysis of the corresponding phase hamiltonian [0 . 
Impurity-induced antiferromagnetism is already exten- 
sively investigated for spin ladders ]l9|-p^ . 

In the next paragraphs we will refine the calculation 
of the norm and the local magnetization by considering 
indeed a discrete wave function. We would like to remove 
the non-orthogonality of our basis states in order to be 
on safer ground for the subsequent reasoning. Moreover, 
certain boundary effects can be captured by this proce- 
dure. 

By Gram-Schmidt orthogonalization wc find the or- 
thonormal basis 



{w,\SI^\w,) = l/3 



\wo) := |0) 
K>o) := (-ir-|(N) + N-l)). 
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The above definitions simplify any norm calculation. The 
magnetization calculation becomes a bit more tedious. 
For odd / = 2j -I- 1 we have 



{W]+i\S2j+i\wj+i) = -1/6 

K-|52,+il^J-) = l/2 forj=0 
Kl^2,+il^j) = 1/3 forj->0 
{wj\S^j+i\w^) = 2-1-^+* for J > i > 



(9a) 
(9b) 
(9c) 
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and the analogous formulae hold for i ^ j. All other 
expectation values are zero. From the relations (|^) one 
finds for a normalized state \v) = X^i^o ^d'^^i) 



{v\S^,+M^\\bA'-]:\b,+i\ 
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(10a) 



(10b) 



{wj\SI^\w,) = -2-1-J+' for < i < j 



From Eqs. (O) we obtain 



(11a) 
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{v\SI,\v) = gifo.f -Re bp-^ l—bo + J2b,T 



(12) 



Note that for I ^ and slowly varying hi one re- 
obtains, of course, the results Eqs.(|) and (0). This con- 
cludes the calculation of the local magnetizations. 



III. KINETIC AND POTENTIAL ENERGY 

In this section we address the important issue of how 
the simple spinous introduced in the previous paragraph 
move and how they are attracted by the chain end. We 
will distinguish two regimes of frustration (i) 0.5 > a > 
ac = 0.241 and (ii) a < a^. 

In regime (i) we know that the spinous have a finite 
mass and that they have a quadratic minimum in their 
dispersion |l^jlj,|l^. Moreover, the strongly frustrated 
chains display spontaneous symmetry breaking of the 
translational symmetry, i.e. spontaneous dimerization 
occurs. The actual value of the gap (and the dimeriza- 
tion), however, is very small up to a = 0.35 |2^j2^ (see 
also Fig.|). 

In regime (ii) the spinon dispersion is linear in the wave 
vector k p^ . Hence we assume uj(k) — vs\k\ where vs is 
the spin wave velocity which is roughly given by 



vs 



-(1 - 1.12a) 



(13) 



according to Ref. ||2^,Q. The multiplication by |fc| in re- 
ciprocal space can be visualized as the consecutive mul- 
tiplications sgn(A:) • k. To multiply by k corresponds to 
— in real space, while a multiplication by sgn(fc) cor- 
responds to the convolution with the principal value of 
i/{'Kx). Then, |fc| in reciprocal space corresponds to the 
convolution 



d 



dx tt{x - y) 



■4'{y)dy (14a) 



= — V 
n 



(x - y)' 



:^{y)dy . (14b) 



Repeating the same for even / = 2j we arrive at 



The symbol V stands for the principal value which has 
to be taken. If one has to discretize the singular operator 
i?kin, the most reliable approach to follow is to go back 
a few steps using uj{k) = vs\sui{ka)\/a in fc-space and 
then transform this to real space yielding 
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(15) 



where a is the lattice constant. 

If we want to treat a chain end, a modification turns 
out to be necessary. Calculating with Eq. ( p^ ) the spinon 
wave function for an undimerized odd size chain one finds 
that the spinon is strongly repelled from the borders. 
This is similar to the problem of having a quantum me- 
chanical particle in a finite box where the wave functions 
are sine-like with nodes at the wall positions. It is due 
to the fact that the kinetic energy cannot be fully satis- 
fied close to a wall since hopping processes through the 
wall are not possible. However, analyzing the numerical 
results (see for instance the uniform susceptibility in 
Fig. 9(b) of Ref. ||l^) one realizes that the probability of 
finding the spinon at a given site is approximately con- 
stant. We interprete here the uniform susceptibility Xi 
as calculated by Laukamp et al. [ p^ as a measure for 
the probability to find the spinon. From this observation 
we conclude that the walls not only truncate the hop- 
ping processes but that they induce also another effect. 
To account for the fact that an approximately constant 
spinon wave function is the ground state wave function 
for the problem without dimerization we modify Eq. ( p^ ) 
to 



{x - y) 



,^{y)dy-^. (16) 



The effect of this modification is most easily understood 
by going back to Eq.(14a) and exchanging the sequence 
of differentiation and convolution. A constant wave func- 
tion yields a (5-function due to the initial jump. The re- 



sulting term is compensated by the 1/x term in Eq.(16). 
We like to stress that the modification is motivated only 
phenomenologically. It would be desirable to have a more 
microscopic justification as well. 

Equivalently, the discrete version of Eq. (|l^) reads 
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if the site counting starts with i = 1 at the border. 
Thus, for chain ends with small frustration we will use 
Eq.(16) as kinetic energy. For bulk problems we will use 
Eq.(14). The corresponding discrete versions are Eq.(p^) 
and Eq.(p^), respectively. 

Thus far, we have discussed only the kinetic energy. 
Considering now the potential energy we first consider 
regime (i) and in particular the Majumdar- Ghosh point 
a = 0.5 for which the ground state is identical to the 
short-range RVB state. As noted previously, e.g. ||l^,|ll|, 
the confining potential is a linear one in this regime. The 
potential due to the dimerization 5Hd with 



is given by 

F(2^ + l) = (5((^|i^D|^)-(l|i^D|l)) 



(19) 



By inspecting Fig.^a) and (c) we find that V{2i + 1) in- 
creases by (3/2) J if the spinon hops once, since one sin- 
glet on the strong bonds is lost on the right side whereas 
one singlet is inserted on the weak bonds on the left side. 
Thus we have 



V{2i + l) =5-Ji 



(20) 



This allows us to propose the following Schrodinger equa- 
tion for the motion of the spinon 

J 9^ 3(5J 
= -'h^dx^'^^''^ + — (21) 

with the restriction a; > 0. The value of m can be found 
from a variational investigation of a single spinon and is 
found to be approximately m = {\ + 7/V65)^^ ~ 0.535 
|l|,0. Rescaling Eq.(|l]) hy x ^ S,y with 

^ = (3m,5/2)-i/3 (22) 

(the site spacing is set to unity) yields 



E^{y) = J 



/ (3V4)^ y/' 



\ 2m 



dy 



^ij{y) + yi^{y)] . (23) 



The linear differential equation ( ^l[ ) with the boundary 
condition '0(0) = is solved by shifted Airy functions 




= -7r2^i(y) + yAi(2/) 



-ZiA\{y + Zi 



dy 
&^ 
dy 



(24a) 



^M{y + Zi)+yki{y + Zi) (24b) 



where the Zj are the zeros of Ai(y) which define by 
ej = —Zi the eigenenergies of the rescaled problem. The 
normalization is given by ■0(2/) = Ai(y-|-z,i)/|Ai'(2/-|-Zi)|. 
Summarizing, using Eq.(p2) we have obtained 



ip{x) 



./Mil: 

Ai(a;/g + Zi) 
Ve|Ai'(2/ + ^.)| 



2\ 1/3 



(25a) 
(25b) 



The form of the wave functions are given in Figj^ where 
^ is made equal to unity. 



2j 



(18) 
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0.0 1.0 2.0 3.0 4.0 5.0 6.0 7.0 8.0 
X 

FIG. 2. Ground state and first excited state wave functions 
ipo{x) and tpi{x) using a linear potential V{x) = x and a 
quadratic kinetic energy //kin = ■ The corresponding eigen 
energies are eo = —20 = 2.338 and ei — —z\ — 4.088. 

A linear potential together with a quadratic kinetic en- 
ergy leads to energies proportional to 5^/^ if 5 measures 
the potential strength. Such a behavior is actually ob- 
served in the dependence of the triplet gap Atrip (^) on 
the dimerization Atrip (5) — Atrip (0) c>c 6"^^^, see Fig.^ and 
p8| . The corresponding calculation is almost equivalent 
to the above one since the triplet state can be viewed as 
two parallel spinous bound together by the same poten- 
tial. Then, the coordinate x corresponds to the relative 
coordinate and the mass m has to be replaced by the 
relative mass fi — m/2 since one has two kinetic energy 
contributions, see e.g. [ p]Jl^ . 

In the light of the success of a simple quantum mechan- 
ical picture for the case of large frustration we come back 
to case (ii) of subcritical frustration. Before entering the 
discussion we emphasize that without a gap there is a 
priori only a less good justification for a one-particle de- 
scription since the length scale of the object carrying the 
S = 1/2 is the same as the spatial extent of the bound 
wave functions. Yet it is interesting to see that the sim- 
ple model works well qualitatively and to a certain extent 
even quantitatively. 

For subcritical frustration, we have two important 
pieces of information. The kinetic energy is linear in 
fc, not quadratic. But the exponent of the (triplet) 
gap growth with dimerization is also 2/3 or close to it 
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j2/3^_ To illustrate this point we 
gap data for four different values of the 
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FIG. 3. Triplet gap as function of the dimerization 5 for 
various values of the frustration a. Symbols are DMRG re- 
sults; the solid lines are simple power law fits A « Ao -I- AiS'^ 
with the parameters {a: Ao, Ai, ly) (0: 0, 1.57, 0.65); (0.241: 
0, 2.05, 0.65); (0.35: 0.033, 2.216, 0.642); (0.5: 0.24, 2.183, 
0.685). The absolute error is 10^* J at most. 

Power law fits with exponents close to 2/3 describe 
fairly well the gap growth in the subcritical and the su- 
percritical frustration regime. The appropriate fit pa- 
rameters, however, depend on the fit interval chosen. It 
is known that logarithmic corrections make it in general 
very difficult to observe the asymptotic behavior. Our es- 
timations for the binding length (see below) provide ad- 
ditional information on how small the dimerization and 
how large the system have to be in order for the asymp- 
totic behavior to be reached. 

Thus far, we keep the information that the confining 
potential cannot be linear in the subcritical frustration 
regime. Revisiting the scaling that mapped Eq.(pl|) to 
the universal form Eq. ( p3| ) shows that the potential must 
grow as a square root V{x) oc ^/x in order to retain the 
exponent 2/3. This seems a very plausible conjecture 
which will be corroborated by a direct calculation of the 
potential. 

In order to evaluate the potential energy between the 
chain end and the spinon we use ideas of Talstra et al. 
P^ , p5[ . These authors have shown that a spinon state 
can be generated to 98% overlap by inserting a single spin 
in a spin chain which is otherwise in its ground state. If 
we take this idea over to our problem of a spinon and a 
chain end we may assume that the systems between the 
chain end and the spinon is in its undimerized ground 
state which is the ground state of a finite piece of chain 
with open boundary condition. We have to calculate the 
expectation value of the dimerization operator Hjy ( [l^ ) 
in this undimerized ground state in order to obtain an 
estimate for the potential. A small refinement is actually 
necessary since we are only interested in the expectation 
value of Hd with respect to the bulk limit. So we set 
V{2i -|- 1) = {H{^)2i where the subscript 2i refers to the 



5 



length of the finite piece of chain and 



{S2j-lS2j — (<S'lS'2)bulk) • 



(26) 



From this equation one sees that the reference to the bulk 
limit yields only a finite offset. 

Note that calculating {H^) with respect to the 
undimerized ground state yielding a potential V{x) oc S 
is in the spirit of degenerate first order perturbation the- 
ory. The effect of the perturbation, here dimerization, 
in the subspace of the unperturbed elementary excita- 
tions is considered. Even though these excitations are 
not really degenerate they are arbitrarily close in energy 
so that none of them can be discarded. 

The calculation of {H{^) is a perfect task for the DMRG 
approach p^j37| |, especially since we deal with open 
boundary conditions. The results are shown in Fig.^. 
We use the infinite system algorithm keeping fOO 
states in each step. The error in the energies due to ba- 
sis truncation is found to be smaller than 10~^. 
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FIG. 4. Expectation value of as explained in the main 
text for three values of a. The solid lines are square root 
fits oc VL - 1.5 with the prefactors 0.505, 0.84, and 1.29 for 
a — 0, 0.241, and 0.35. The offset 1.5 is chosen for the im- 
provement of the fits. The dashed lines are power law fits 
oc (L — 1.5)'' with the exponents and prefactors 0.544, 0.431 
for a = 0, 0.630, 0.545 for a = 0.241, and 0.804, 0.450 for 
a = 0.35. Inset: Square of the same data for longer chains. 

One clearly observes a sublinear increase of {H{^) with 
the chain length. Without frustration (a = 0) a square 
root power law fits the data perfectly if a small offset on 
the X-axis is taken into account. For larger frustration 



larger exponents yield better fits. The inset, however, 
shows that for longer chain lengths (L > 50), the square 
root behavior is recovered for a < ac which can be seen 
from the linear behavior of the squared values. The frus- 
tration value a — 0.35 is included to show that for this 
value relatively close to the critical one a linear behavior 
(H^) cx L cannot yet be seen. 

Another evidence that for a < ac the expectation value 
(TJq) rises proportional to a/L stems from the dimension 
1/2 of the dimerization operator By integration 

/ dx/^/x cx ^/L follows the conjectured behaviour. 

So our conjecture of a square root confining potential 
for lower frustrations is corroborated by direct calcula- 
tions. We are now in the position to perform an anal- 
ysis as in Eqs.(^ -p5|). Starting from the generalized 
Schrodinger equation with the border effect corrected ki- 
netic part H^.^^ Eq.® 

(27) 



we rescale hy x — with 



\saj 



(28) 



where vg is the spin wave velocity and A some con- 
stant which can be deduced from fits to data as in Figj^. 
Eq.(p7|) becomes 



from which the energies are found using 



(29) 



(30) 



The resulting two first wave functions are shown in 
Fig.^. The dimensionless energies are eo = 1.049 and 
ei = 2.040. Equation ( ^9| ) is solved numerically by 
equidistant discretization and fast Fourier transform with 
up to 2^^ points. Due to the singular behavior of Eq.(|29|) 
the precision is only of about 10"'^. But this is sufficient 
for the intuitive description that is being proposed here. 
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FIG. 5. Ground state and first excited state wave functions 
^po{x) and ^pi{x) using a square root potential V{x) = ^/x and 
linear kinetic energy -ffkin = \k\ plus border corrections (see 
Eq.(p^)). The corresponding energies are eo — 1.049 and 
ei = 2.040. 

Note that due to the border correction Eq.(p^ the 
wave functions do not vanish at the chain end, i.e. at 
X = 0. Furthermore, the wave functions decay rapidly 
as X oo, but not following an exponential form, e.g. 
1^,^, but rather a power law il>{x) oc This is due 

to the highly non-local character of the kinetic energy 
Eq. ( p^ . Thereby, we do not want to claim that the decay 
for large x follows indeed a power law. The picture we 
are presenting does not contain the generation of further 
spinon-antispinon pairs since it is a quantum mechanical, 
single particle scenario. It is the pair creation possible 
for energies larger than the gap Atrip which eventually 
leads to the exponential decay of the spinon probability. 
Our picture works best at intermediate distances as we 
will show in comparison to the DMRG data in the next 
section. 



IV. RESULTS: CHAIN ENDS 

Now we are in the position to compare results of bound 
spinon calculations directly against numerical results. 



A_ 

N 

V 



V 



A_ 

N 

V 



0.25 
0.15 
0.05 
-0.05 
-0.15 
0.15 
0.05 
-0.05 
-0.15 
0.15 
0.05 
-0.05 
-0.15 



□ 


1 ■ 1 ■ 1 ' 

a=0 

8=0.048 






1 




1 1 1 1 ■ 1 ■ 1 1 

a=0.241 
8=0.012 

I I . 


' 1 ' 

'^"^^^ 


1 ■ 1 ■ 1 ' 

a=0.35 
8=0.012 


1 1 i_ 





10 



20 30 

site i 



40 50 



FIG. 6. Local magnetizations at dimerized chain ends with 
weak bonds for various values of a and S. St^t is set to 1/2. 
Squares: odd sites; circles: even sites; filled symbol: DMRG; 
open symbols: discrete bound spinon calculation; dashed lines 
continuous bound spinon calculation. 



The latter are obtained using the DMRG technique 
[p6| , p7t for open boundary conditions with between m = 
24 to 32 states kept in the iterations. The results are 
shown in Fig.^. Quantum mechanical bound spinon cal- 
culations are carried out in two ways. One is the contin- 
uum one relying on Eqs.(^,^ 16), the input for the spin 
wave velocity Eq. (|l^) and the potential parameters for a 
square root potential (see caption of Fig.H). This means 
that one takes the result for |V'(a;)p in Figl^ and rescales 
iV'Ca:^)!^ with ^ from Eq. (^) . The relevant 

values are ^ = 16.1, 23.5, and 15.6 for the uppermost to 
the lowermost panel in Figj6[ The results are shown as 
dashed curves. 

The other one is a discrete treatment relying on Eqs. 
( |lO| , p^ , p7| ) . The results are depicted with open symbols. 
The values of the potential are those read off FigJJ with 
V{2i + 1) = {H{^{L — 2i)). For simplicity the very good 
power law fits (dashed lines in Fig.^ are used for the 
numerics. 

The agreement between the DMRG results and the 
bound spinon model is very good, especially in view of 
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the simphcity of the model. The shape of the curves are 
very well described by the bound spinon results. In par- 
ticular, the correct binding length is predicted. This is 
in particular interesting since we know that the quantum 
mechanical model does not display an exponential de- 
cay. This means that the exponential tail matters only at 
larger lengths. For intermediate lengths, where the main 
weight is found, the bound spinon calculations work very 
well. The main feature that is included by the disrete cal- 
culation is the decrease of the amplitude for the first four 
to five sites. It results from the properties of Eqs.(|0 12) 
and can thus be attributed to border effects which were 
neglected in Eqs.(^,0) but properly taken into account in 
Eqs.(0|2|). It is not due to the dynamics of the problem. 
The very good agreement of the open-symbol curves and 
the long-dashed curves for larger distances underlines the 
applicability of the approximation (^. 

The overall amplitude is larger than the one predicted 
by the bound spinon model. Since the basis we used 
included only the shortest range singlets it is not sur- 
prising that the true antiferromagnetic correlations are 
in fact larger. Thus we view the bound spinon results in 
this respect as a lower bound for the antiferromagnetic 
correlations. But it is interesting that even for the ampli- 
tudes the agreement is reasonable for correlation lengths 
^ around 15. For larger values of ^ (middle panel in FigJ^) 
the agreement becomes not as good. 

Another point that can be addressed easily within the 
bound spinon model are excitation energies, namely the 
energy difference Ei — Eq (^) between the ground state 
-00 and the first excitetd state ^i, see Fig.| The con- 
tinuum (discrete) calculation yields 0.97 (0.10), 0.048 
(0.055), and 0.061 (0.073) for the three cases in Fig.| 
in descending order. It is not surprising that the contin- 
uum calculation provides lower values since it assumes a 
potential with smaller exponent, see Fig.^, so that the 
distances between the energy levels become smaller. A 
smaller exponent implies that the potential increase for 
smaller arguments is larger than for larger arguments. 
Thus the smaller of two consecutive eigenenergies is lifted 
with respect to the larger one since the wave function be- 
longing to the smaller eigenenergy is more localized. This 
effect explains also why the deviation is the largest for 
a = 0.35 where the deviation of the optimum fit expo- 
nent to 1/2 is largest. 

It was mentioned before that the occurrence of the 
triplet gap in dimerized systems can be viewed also as 
a binding phenomenon, see e.g. [ pT|Jl3| . In this case, two 
spinous bind. However, since only the relative coordinate 
matters, this problem is almost equivalent to the binding 
of a single spinon to a chain end. The main difference is 
that the kinetic energy is doubled since both interaction 
partners move. So it is understandable that the ratio 
{El — i?o)/ Atrip is fairly constant. The triplet gaps Atrip 
as read off Fig!| are 0.225, 0.116, and 0.162. The ratios 
of the continuum (discrete) values are 0.43 (0.46), 0.42 



(0.47), and 0.37 (0.45). Similar to the result found pre- 
viously (0.6) ||l^ the transition from the ground state to 
the first excited one is at about half the triplet gap. This 
implies that it should be observable in scattering exper- 
iments as a rather sharp feature within the gap. Indeed, 
Raman scattering results have revealed such a feature in 
Zn-doped CuGeOs even though the energy is higher 
than the one-dimensional model predicts. 

It is clear that at energies Eq + Atrip pair production 
becomes possible as it is known in QCD. The gedanken 
experiment of separating quarks in space in spite of the 
confinement results in the production of new quark pairs. 
In our situation any excited bound state can decay by 
producing a pair of spinous at low values of momentum 
if its energy fulfills Ei > Eq + Atrip. This mechanism 
gives rise to a continuum which sets in at ii^o+Atrip. Thus 
higher energy bound states are no longer true eigenstates 
of the problem but have a finite lifetime since they might 
decay. The discrete calculation actually shows that E4 
and higher are no longer stable. This means that be- 
sides El, also E2 and E3 should exist as distinct sharp 
modes. Experimentally there is so far no indication for 
these higher modes. But it should be borne in mind that 
any deviation from pure d — I behavior tends to decrease 
the binding energies. Bound states are shifted closer to 
the continuum or even disappear in the continuum due 
to higher dimensional effects, see also ]3^Jl^. 
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FIG. 7. Same as Fig.^ but for a = 0.5 and 5 = 0.012. 

In Fig.^ we show the same results as in Fig.^ but for 
the Majumdar-Ghosh point a = 0.5. They are based 
on Eqs.(pTU^|7|) in the continuum case. In the discrete 
case we used Eqs.( p0| , p^ , p^ . The discrete kinetic energy 
used is nearest neighbor hopping such that the quadratic 
minimum is the same as for Eq.(^). 

One should expect an improved agreement since the 
kind of basis we are working with is particularly suited 
for the Majumdar-Ghosh point |jl^,|l^. The agreement, 
however, is not very good. It becomes decisively better if 
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the bound spinon results are shifted by about three lat- 
tice spacing. Probably dynamic border effects not taken 
into consideration by the basic assumptions concerning 
the kinetic and the potential energy (focused on the low 
energy behavior) are important. The binding length ^ for 
the Majumdar- Ghosh model according to Eg. ([2^) is only 
4.7. Thus a low-energy approach might be insufficient for 
such a small binding length. 

Interestingly, Eggert and Affleck found in a compari- 
son of continuum field theoretical results with numerical 
data at ^ = also that the agreement is much better if 
the continuum results are shifted by two sites towards 
the chain end |39|. An understanding of either of these 
observations would help to understand the other. 



V. RESULTS: KINK DEFECTS 

Spin-Peierls systems are characterized by a coupling of 
magnetic (quasi) one-dimensional spin degrees of freedom 
to the lattice degrees of freedom, i.e. vibrating distortions 
(for reviews see pO| , ^ ). At low temperature and zero 
magnetic field the system is in the gapped dimerized (D) 
phase where the coupling strength Ji of the spin chain 
alternates from site to site. As the magnetic field is in- 
creased the gap becomes smaller and at a critical field 
the phase changes to an incommensurably modulated (I) 
phase. This I phase can be viewed as a soliton lattice 
where equally spaced tanh-like zeros of the distortions 
5i occur. In the vicinity of each zero of the distortion 
a spinon S* = 1/2 is localized, see and references 
therein. 

In Fig.^ the soliton lattice is illustrated. These results 
are found from an adiabatic calculation where the distor- 
tions are treated statically. Most of the investigations of 
the I phase assume static distortions. The Hamiltonian 
H from which Figjs] is computed by minimization of the 
ground state energy with respect to {5i\ reads 



-Bclast (31a) 

'aS^-S,+2) (31b) 

(31c) 
(31d) 

(31e) 



where a denotes the relative frustration and Sz is the 
z component of the total spin of the L-site chain. The 
last two terms in Eq. (|3^) are the Zeeman energy and the 
elastic energy associated to the lattice distortion. A site 
independent spring constant K is used. 
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FIG. 8. Two solitons in a hundred site chain, i.e. 
m = 1/100, with periodic boundary conditions are studied 
(parameters a — 0.35, K — 8.7J). In the upper panel, re- 
sults from the self-consistent determined distortions Si (cir- 
cles) which minimize {H) in (|^) are shown. They are de- 
picted between the sites since they belong to bonds. The 
straight lines are fits to the linear rising distortions with slope 
s = 0.0045. The lower panel displays the resulting local mag- 
netizations (circles, DMRG calculation) and the bound spinon 
result (solid lines). For details see main text. 

There are regions around site, e.g., 50 where the dis- 
tortion alternates as in the D phase. The corresponding 
local magnetization is essentially zero. Close to the linear 
zeros of the distortions at about sites 25 and 75 one ob- 
serves a strongly alternating local magnetization which 
adds up to 1/2: T,Zoi^i) = V2- The distribution of lo- 
cal magnetizations is observable by NMR and gives 
interesting information on the I phase and on the validity 
of the adiabatic treatment Q]. An interesting observa- 
tion is that the ratio between the parallel magnetizations 
to the antiparallel magnetization at (roughly) the same 
sites is 7:4 as implied by Eqs.(||j^). 

In this paper we are not interested in the self-consistent 
minimization of (^ij) but we take the {Si} as given and 
use the quantum mechanical picture developed above to 
show that the localized alternating magnetization in the 
vicinity of a zero of the distortion results from a bound 
spinon. The zero of the distortion is viewed as a kink 
defect. Since the 6i on the odd (even) bonds change sign 
at the zero there is no simple short-range singlet pattern 
which ensures singlets at the stronger bonds without any 
free spin. Insofar the situation is comparable to the one 
at chain ends with weak bonds at the ends. The spinon 
bound to a kink defect is even simpler than the spinon 
bound to a chain end since no particular boundary effects 
on the kinetic energy need to be taken into account. Thus 
we use the continuum description based on Eqs.(^,0,^. 

The derivation of the attractive potential requires an 
additional assumption in order to treat also site depen- 
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dent distortions. We find it natural to assume that the 
chain end potential V{x) as determined from Fig. ^results 
from expectation values (SiSi+i) — {SiSi+i)huik which de- 
cay like a power law with increasing distance from the 
ends 



V{x) = Sax'' 
ah 



{Sy''-^ + 5{x - vf-^)dy , 



(32) 



where we took into account that the chain end and the 
position of spinon constitute borders. The effect of these 
borders is treated approximatively as additive, i.e. as 
independent. A general power law with exponent b and 
prefactor a is dealt with. The advantage of Eq.(^) is 
that the distortion can be made site dependent d —^ 6{x). 
In the vicinity of the kink defect, d{x) — sx is a good 
description where s is an appropriately determined slope, 
see upper panel in Fig.|[ A short calculation then leads 
to 



V{x] 



kink 



abs 
as . 



x\ 



b+1 



iyy'>-'+y{\x\-y)''-')dy 



(33) 




-4.0 -3.0 -2.0 -1.0 0.0 1.0 2.0 3.0 4.0 
X 

FIG. 9. Ground state and first excited state wave functions 
ipo{x) and tpi{x) using the integrated square root potential 
V(x) = and the linear kinetic energy Hkin = \k\ (see 

Eq.(|l4)). The corresponding energies are eo ~ 1.051 and 
ei = 2.310. 

Thus we obtain again a power law confining potential, 
the exponent of which is increased by one compared with 
the site-independent distortion. Moreover, the potential 
does not depend on whether the spinon moves to the 
right or to the left of the zero. 

In Fig.| a generic result for the localized state at a 
kink defect is shown. All parameters are set to unity. 
But rescaling permits, as before, to obtain the general 
relations 



/2vs_ 
\ as 



l/{2+b) 



E,, 



,1+b' 



as X 1/(2+6) 



(34a) 
(34b) 



The solid lines in the lower panel in Fig.^ depict the 
convincingly good agreement between the bound spinon 
picture and the DMRG calculation. The value s = 0.0045 
and the parameters of the dashed line fit to the a = .35 
data in FigJJ is used for the binding length <^ in (34a). 



VI. SUMMARY 

In this work we have developed a quantum mechanical 
picture in the spirit of first order degenerate perturbation 
theory which describes the confinement of spinous due to 
dimerization. Without dimerization there is no confine- 
ment independently of the value of the frustration a. The 
picture works below and above the critical frustration of 
the undimerized chain ac = 0.241, even though the jus- 
tification below ac is less good. The main ingredients 
of the calculation are the known facts concerning the ki- 
netic energy of the spinous which is linear below ac and 
quadratic (plus a mass term) above ac- From the dimer 
expectation values close to the borders of finite chains 
we deduced a potential V{x) which is proportional to 
the dimerization S. It was shown that this potential can 
be nicely fitted for intermediate distances (L w 50) by 
sublinear power laws below and just above ac = 0.241. 
For larger distances (L — > oo) our results indicate that 
V{x) increases like a square root for a < ac and linearly 
for a > ac- To our knowledge, the fact that the confining 
potential below ac = 0.241 has to be sublinear has not 
been reported before in the literature. 

The quantum mechanical picture describes very well 
the local magnetizations close to chain ends with weak 
bonds for intermediate ^ (« 20). A very good agree- 
ment is also found for the spinon bound to a kink de- 
fect. The sum of the moduli of the local magnetizations 
is (at least) enhanced by a factor 11/3 over the value 
for a single spin. For low ^ values the continuum de- 
scription becomes less reliable. For large values of ^ the 
form of the site dependence of the local magnetization 
is still captured by our quantum mechanical picture. In 
particular, binding lengths can be easily estimated using 
Eqs.(22,28 34a). The overall amplitude of the staggered 
component, however, is in fact larger for 5 — > at T = 
11 ■ 

Besides local magnetizations the eigenenergies of the 
bound states can be estimated easily within our ap- 
proach. We find that in truly one-dimensional systems 
transitions to the first three excited states at chain ends 
should be possible. Any degree of higher dimensionality, 
however, will reduce this number. 
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The formation of gapful triplets can equally be viewed 
as binding of two spinons Q. The essential differences 
to the case of chain ends are that the site variable denotes 
the relative coordinate between the spinons and that the 
kinetic energy is doubled since both interaction partners 
move. 

The quantum mechanical picture is scale invariant. 
This means that for power law potentials the wave func- 
tions and energies can all be scaled to one universal case, 
see Figs.^ and ||. Thus the quantum mechanical de- 
scription provides a pedestrian approach to perturbations 

with anomalous dimensions and exponents. The anoma- 
lous exponents result naturally from the scaling behavior 
of the potential and the kinetic energy under changes of 
the length scale. A dimensional analysis implies already 
that for V{r) = ar^ and i?kin = c|fc|'' the characteristic ^ ' 



lengths scale like f oc (c/a)^/'^''"'"''^ and the characteristic 
excitation energies like E oc {a'^c^Y^'^^^'^\ 

We like to point out that the basis states of our ap- 
proach can be refined to take longer range antiferromag- 
netic correlations better into account. To this end longer 
range singlet pairs must be included and superposed. 
Dispersions and potentials can then be found by con- 
cepts similar to those proposed by Sutherland for 
short range singlets in d = 2. 




FIG. 10. Sketch of the singlet distribution on a doped 
ladder with strong bonds (thick solid lines) on the rungs. The 
open eyelets represent singlets. The circle indicates the defect 
(nonmagnetic dopant), the arrow the generated spinon after 
three hops. The misaligned diagonal singlets generate the 
confining potential. 

As an outlook to higher dimensions we present in 
Fig.|l^ the leading sketch of the singlet distribution on 
a ladder with strong bonds on the rungs. Ladders are 
objects of ample investigations. The alternating magne- 
17,^^ and enhanced antiferromagnetic cor- 
indicate the similarity to open chains. 



tizations 
relations 



In Fig.hdthe situation is depicted of the spinon (arrow) 
generated by the insertion of the defect (circle) after three 
hops on the upper leg. One clearly sees that the motion of 
such a spinon is also confined by a monotonic increasing 
potential. An increasing distance of the spinon from its 
original rung induces an increasing number of misplaced 
singlets. This leads to an increase in energy. From Fig.p^ 
it is clear that the basic concepts of spinon motion and 
confining potential can be extended also to other RVB- 
type spin systems. The results discussed in this paper 
are not restricted to chains only. 
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